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ABSTRACT 

The present study is an attempt to investigate the long term variations in coronal 
rotation by analyzing the time series of the solar radio emission data at 2.8 GHz fre- 
quency for the period 1947 - 2009. Here, daily adjusted radio flux (known as Penticton 
flux) data are used. The autocorrelation analysis shows that the rotation period varies 
between 19.0 to 29.5 sidereal days (mean sidereal rotation period is 24.3 days). This 
variation in the coronal rotation period shows evidence of two components in the vari- 
ation; (1) 22-years component which may be related to the solar magnetic field reversal 
cycle or Hale's cycle, and (3) a component which is irregular in nature, but dominates 
over the other components. The crosscorrelation analysis between the annual average 
sunspots number and the coronal rotation period also shows evidence of its correlation 
with the 22-years Hale's cycle. The 22-years component is found to be almost in phase 
with the corresponding periodicities in the variation of the sunspots number. 
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1 INTRODUCTION 

Coronal rotation can be observed through various solar trac- 
ers at different frequencies, like coronal green line (Fe XIV 
emission line at 530.3 nm ), white light, He I line (at 1083 
nm), soft X-rays, UV rays, radio waves. The coronal green 
line has been used to measure t he rotation r ate o f the so- 
lar co ro na at higher la t itude s by Waldemierl (|l950l); Trellis 



19891) : lRvbaklll 19941) -iBadalvan et all ll2006al. 



1957 1) : ICooper et al.1 (Il962l); ISvkoral (Il97ll ): ISime et al 



- and oth ers. 

The results of l Waldemierl l|l950l ) and lCooper et al.1 j 19621 ) in- 



dicate a faster rate of rotation as compared to the rate of 
rotation of the sunspots, suggesting a much lower differen- 
tial rot a tion r ate in the corona. In his work on green corona 
ISvkoral (|l97ll ) found that the Sun shows little or no dif- 
ferential rotation for six latitudinal zones ±7.5, ±27.5 and 
±47.5. For lo w latitudes th e rotation period was near to 
that found by iTrellisl j 19571 ). The green (Fe XIV at 530.3 
nm) emission line for the period 1973 - 2000 and red (Fe X 
at 637.4 nm ) emission line for th e period 1984 - 2000 were 
analyzed bv Ultrockl (|l997l . 120031 ). It was reported that the 
corona, at green and red emission l ines, shows more rigid 
rotation than does the photosphere. ISime et alj (1989) also 
concluded, after analyzing the Sacramento Peak Observa- 
tory (SPO) data observed between 1973 and 1985, that the 
Fe XIV corona rotates more rigidly than do features in the 
photosphere or chromosphere. The synodic period obtained 
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bv lRvbakl | 1994) for the period 1964 - 89 again confirmed 
the diff erential rotation of the green corona. [Badalvan et al.l 
( 2006a. b) carried out a comprehensive analysis using a long 
database (1939 - 2001) on the brightness of the coronal green 
line. The results support previous conclusions that the dif- 
ferential rotation in the corona is less pronounced than in 
photospheric tracers. 

lHansen et al.l il 19691 ) used the K-coronometer for coro- 
nal rotation measurement at different latitudes, for heights 
ranging from 1.125 to 2 Rq. The rotation found at the equa- 
tor is in good agreement with the sunspot 's rotation results 
and shows less variation with latitude at higher latitudes in 
comparison to the rotation of the chromosphere. A detailed 
study of the white light corona, from 1.1 _Rq to 30 -Rq, was 
done with the LASCO coronagraphs on board the SoHO 
spacecraft. It was concluded that the rotation of the corona 
displayed a radially rigid rotati on of 27.5 days sy nodic pe- 
riod from 2.5 Rq to > 15 Rq l|Lewis et al.lll999l ). 

The He I 1083 nm maps, from the National Solar Ob- 
servatory (NSO), have been used to determine the rotation. 
It is found, both from observations and from magnetic ex- 
trapolation methods, that the corona becomes more rigid 
with height. By considering coronal holes as tracers (from 
an atlas of coronal holes mapped in He I 1083 nm data) 
of the differential rotation llnslev et al.l (|l995l ) demonstrated 
that the mid-latitude corona rotates more rigidly than the 
photosphere, but still exhibits significant differential rota- 
tion, with an equatorial rate of 13.30 ± 0.04 deg/day, and 
at 45°latitude a rate of 12.57 ± 0.13 deg/day. An analy- 
sis of the rotation of coronal holes spanning 18 years (from 
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Figure 1. The annual mean of solar radio flux at 2.8 GHz is compared with the annual mean of sunspot numbers for the period 
1947-2009. 



1973 to 1991) wa s done based on data from the Catalogu e 
of Coronal Holes (|Navarro-Peralta fc Sanchez-Ibarrall 19941 ). 
Isolated coronal holes showed a typical differential rota- 
tion, but polar coronal holes extensions displayed two differ- 
ent types of behavior: a rotation rate below approximately 
40° ± 5° of hcliographic latitude, increasing to the equator, 
and a rotation rate above the same heliographic latitude but 
increasing to the poles. 

Coronal holes, as observed from the Skylab and Yohkoh 
spacecrafts, have also been used to determine the rotation 
rate of the outer corona. Soft X-ray (SXR) observations 
of an elongated coronal hole, shows the almost rigid rota- 
tion of the coronal hole (|Timothv et al ] |l975l : iKozuka et al] 
1994). The solar full disc (SFD) images, obtained by soft 
X-ray telescope (SXT) of Yohkoh space observatory, were 
used by different scientific groups to study the rotation 
rate of the corona. I Weber et al. I (ll999T l; IWeber fc Sturrock I 
(2002) and IChandra et al. I l|201of concluded, after analyz- 
ing SXT data by different methods, that the rotation pro- 
file of the corona across the latitude is shallow er than the 
rotati on profile of its lower atmospheric levels. iKarivappal 
(2008) tracked the X-ray bright points (XBPs) on SFD im- 
ages observed through the SXT and XRT (X-ray telescope) 
on board the Yo hkoh and Hinode spacecrafts, respectively. 
IKarivappal (|2008T ) found, contrary to all expectations, that 
the corona rotates differentially with respect to latitude, as 
in the case of the photosphere and the chromosphere. 

iKarachik et al. (2006) analyzed the coronal bright 
points (CBPs) on SFD filtergrams observed through 
SoHO /EIT (Fe XII line at 19.5 nm ) and reported that 
the rotation of CBPs closely follows the latitudinal rota- 
tion profile of the photospheric magnetic field. It was also 
shown that coronal features at di fferent heights in the coron a 
exhibits diff e rent r otati on rates. Braisa et alj (|2002l . 2004): 
iMulec et all (120071 ) and lBraisa et al] ^OOsf idetermined the 
solar differential rotation by tracing coronal bright points 
on SFD filtergrams observed through SoHO /EYI (Fe XV 



line at 28.4 n m ). F or the declining phase of solar cycle 23, 
IZaatri et af] (2009) compared differential rotation of sub- 
photospheric layers derived from GONG++ dopplergrams 
with the small bright coronal structures (SBCS) observed 
through SoHO /EIT . It is found at the equator that the 
SBCS rotate faster than the upper subphotospheric layer 
(3Mm) by about 0.5 deg/day. The latitude gradient of the 
rotation rate of the SBCS and the subphotospheric layers 
were found in agreement. 

iKane et al] (|200ll ); IVats et alj (120011 ) used the disk in- 
tegrated solar radio fluxes for the period 1997 to 1999 at 
different fre quencies to study the coronal rotation at differ- 
ent heights. IVats et ahl (|200ll ) shown that coronal rotation 
depends on the heights in the corona. The autocorrelation 
analysis of radio SFD images at 17 GHz, for the period 1999 
- 2001, also ascertained that the solar corona rotates less dif- 
ferentially than does the photosphere and the chromosphere 
jChandraet~ai]|2009l). 



Rvbakl l| 19941 ) reported that the differential rotation of 



the g reen corona shows no clear cyclic variation of the rota- 
tion. IKarivappal (|2008l ). who tracks the XBPs on the SFD 
images of Yohkoh /SXT and i/inode/XRT, also found that 
the rate of coronal rotation is independent of the phases 
of the solar activity cycle. But, from the GONG (Global 
Oscillation N etwork Group) an d MDI (Michelson Doppler 
Imager) data. Ijain et alj (|200ll ). have shown that the equa- 
torial rotation rates show var iation with the solar activity 
cycle. iBadalvan et alj (|2006al lbT) also reported that there is 
periodic variation in the rate of differential rotation of green 
corona which is related to the solar activity cycle, a correla- 
tion usually found in photospheric tracers (jHathawav et al] 
Il990l ). With the SF D filtergrams obtaine d from SoHO /EIT 
(Fe XV at 28.4 nm). lBraisa et al] \2004 ) exhibited that the 
differential rotation pr ofile corresponds ro ughl y to the rota- 
tion o f sunspot groups. IMulec et all l|2007l ) and lBraisa et al] 
(2008) found, with the same data of different epoch, that the 
equatorial rotation velocity and the gradient of differential 
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rotation show similar values during periods of low and high 
acti vity. 

iMouradian et alJ l|2002h analyzed the 2.8 GHz radio 
emission flux with the maximum entropy method and 
showed that the 2.8 GHz radio e mission r otatio n varies ac- 
cording to the activity level. But Mehtal (|2005T l. after ana- 
lyzing the radio emission data could not find any systematic 
relationship between th e coronal rotat i on period a nd the 
phase of the solar cycle. I Chandra" et al l (|2009l . |201Ch calcu- 
lated the equatorial rotation rate for each year separately 
(using radio images at 17 GHz for the years 1999-2001 and 
SXT data for the years 1992-2001, respectively) and com- 
pared it with the annual sunspot numbers. The comparison 
shows that the equatorial rotation rate seems largely to be 
a fu nction of the pha ses of the solar cycle. 

IJavaraiahl l|200Ch used the GPR (Greenwich Photoheli- 
ographic Results) data on sunspot groups compiled during 
1879 - 1975 to study the variations of the differential rotation 
coefficients A and B during the odd numbered solar cycles 
(ONSCs) and during the even numbered solar cycles (EN- 
SCs). The parameters A and B are measures of the equa- 
torial rotation rate and latitude gradient of rotation rate, 
respectively. It is seen that the variation in A is significant 
only in the ONSCs (Waldmeier cycle number 13, 15, 17 and 
19). Whereas, the variation in B is quite significant in both 
ONSCs and ENSCs (12, 14, 16, 18 and 20). There exists a 
good anticorrelation between the mean variation of B during 
ONSCs and ENSCs, suggesting the existence of a 22-years 
periodicity in B. 

Different methods and different tracers have been em- 
ployed to determine the coronal rotation and its dependence 
on the phases of the solar cycle, but the estimates were not 
found to be in agreement with each other. Using radio mea- 
surements to estimate solar rotation is not very new, but 
since the last couple of decades many scientific groups have 
employed radio flux da ta at various frequencies to ascer- 
tain the rotation perio d (|Vats et al II 1998 a; Ka ne et al ] |200ll ; 
IMouradian et af]|2002i'l. its differential ity as a function of al - 
titude (|Vats et alj|2001f ) and latitude jChandra et alj|2009h 
and i ts relationship with the phases of the solar cycle ( Mehtal 
2005). The present study is an attempt to estimate the coro- 
nal rotation rate and its relationship with the solar magnetic 
cycle by analyzing the time series of the radio emission data 
at 2.8 GHz frequency (wavelength 10.7 cm) for the period 
1947 - 2009. The data covers almost six solar activity cycles 
and is currently one of the best and longest data of solar in- 
dices we have after the sunspot number. This radio emission 
at 2.8 GHz frequency originates de ep in the solar co rona and 
hence represents coronal rotation (Vats et al. 200lT). 

The next section discusses the observation and impor- 
tance of radio flux data at 2.8 GHz used in the present anal- 
ysis. Its subsection briefly describes the autocorrelation and 
crosscorrelation analysis technique, primarily used to deter- 
mine the rotation period and its correlation with the solar 
activity, respectively. The paper ends with a discussion of 
the results obtained. 



2 DATA 

The data used in the present work comprises disc inte- 
grated flux measurements of radio waves at 2.8 GHz. It is 



expressed in solar flux units (1 SFU =10~ 22 Watts/meter 
square-Hertz). It is almost thermal in origin because emis- 
sion is mainly due to therma l free-free (bremsstrahlung) ra- 
diation (|Tapping et al.f l990l . It is directly related to the to- 
tal amount of magnetic flux. The data measurements have 
been continued through the Solar Radio Monitoring Pro- 
gramme, a service which is operated by the National Re- 
search Council of Canada. The solar flux density at 2.8 GHz 
is measured using two radio telescopes (called flux moni- 
tors). The two instruments record the strength of the solar 
radio emission at 2.8 GHz frequency each day for as long 
as the Sun is above the horizon. Between 1946 and 1990, 
the measurements were made in the Ottawa area at the 
Algonquin Radio Observatory. In 1990, following the clo- 
sure of that observatory, the system was relocated to the 
Dominion Radio Astrophysical Observatory, near Penticton, 
British Columbia. It now forms a consistent, uninterrupted 
data base covering more than 63 years. The observations, 
which started in late 1946, are still continuing. The conti- 
nuity and consistency of calibration of these measurements 
of the 2.8 GHz radio flux have helped to establish it as an 
internationally accepted index of solar activity. 

These data are readily available through the National 
Geophysical Data Center (NGDC), NOAA, USA. The data 
available at the site of NGDC are tabulated in two forms; 
the observed flux and the adjusted flux. The observed flux 
data are the actual measured values and are affected by the 
changing distance between the Earth and Sun, whereas the 
adjusted flux data are scaled to the standard distance of 
1 AU. Therefore, the adjusted flux data are considered to 
be more appropriate for the study of the behavior of the 
solar corona. The annual mean of the daily measurements 
of the radio emission can be represented in the form of a 
time series (shown by a solid line in Figure [1]). This figure 
shows the yearly variation of sunspots number as well as the 
radio emission at 2.8 GHz. 

ICovingtonl <|l969T l made a comparison over more than 
a solar activity cycle and showed that there is indeed a 
linear correlation between the 2.8 GHz solar flux and the 
total photospheric magnetic flux in the active region. Vari- 
ous manifestations of solar activity are driven by the total 
amount of magnetic flux emerging through the photosphere 
into the chromosphere and the corona, and its spatial and 
temporal distribution. Besides this, the 2.8 GHz solar flux 
correlates well with some other indices of solar activity such 
as sunspot number and total sunspot area. The radio mea- 
surements have an advantage over other indices in that these 
measurements are completely objective, and can be made 
under almost any weather conditions. Since they are closely 
correlated with the magnetic activity which modulates the 
Sun's energy output with solar irradia nce, they correla te 
closely with other solar activity indices (| Covington! fl969V 

The fractal analysis of the solar radio emissions in the 
frequency range from 245 MHz to 15 GHz has established 
that the 2.8 GHz is the most appr opriate frequency for the 
study of the solar coronal rotation (Vats et al. 1998b; Mehtal 
2005). It is so because the fractal dimension is found to be 
the least at a frequ ency near 3 GHz. T herefore, according to 
IVats et all (|l998rj ) and lMehtal (|2005l ). the time series of ra- 
dio emission data available at 2.8 GHz frequency must have 
a very strong modulation due to the solar coronal rotation, 
in comparison to other frequencies. Thus, they used radio 
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emission at 2.8 GHz frequency to determine the coronal ro- 
tation period. We, therefore, used this radio frequency to 
study long term variation in the coronal rotation period, so 
obtained. 

From the knowledge of the yearly mean sunspot num- 
bers (International Sunspot Number, credit: SIDC, 1947- 
2009) on the photosphere, the years 1954, 1964, 1976, 1986, 
1996 and 2009 are the years of solar minima in the solar 
cycles numbered as 18, 19, 20, 21, 22 and 23, respectively 
(Figure [T]). In fact the year 2006 would have been a solar 
minima and the number of sunspots would have increased 
after that, but surprisingly in cycle 23 the minima has be- 
come a much extended one. The sunspots number showed a 
decreasing trend up to the end of the year 2009. It is only 
now in 2010 that solar activity shows signs of increase. We 
can see that the radio flux at 2.8 GHz frequency clearly ex- 
hibits the 11-years cycle performed generally by the sunspots 
number on the photosphere. It is well established evidence 
of a good correlation between the sunspots and the radio 
emission. 



3 METHODS 

There are many statistical methods to analyze such periodic 
time series. In the present work, the autocorrelation and 
crosscorrelation method are employed for the analysis of the 
time series of radio flux at 2.8 GHz. 

3.1 Autocorrelation Analysis 

A time series is defined here as a sequential collection of data 
observations indexed over time. The series of autocorrelation 
coefficients represents the correlation between the successive 
observations of a single time series. The daily measurement 
of radio flux at 2.8 GHz frequency can be represented in the 
form of a time series, (xo, xi,X2, X3, x n -i). 

Now we can form (n — /) pairs of observations from 
this series, where n is the length of the data series (365/366 
days in our analysis) and I is referred to as the lag. They 
are (x ,xi), (xi, xi+i), (x 2 ,X2+i), (x n -i-i,x n -i). 

By regarding the first observation on each pair as one 
variable and the second observation on each pair as the sec- 
ond variable the autocorrelation coefficient P x (l) as a func- 
tion of the lag I is given by 

p.(q = p.h) = E^fr* ~ w ^ Xk+ ' (i) 

where x is the mean of the daily observations. 

The daily solar flux at 2.8 GHz is converted into time se- 
ries of one year each (total 63 data sets) for the period 1947- 
2009, each time series beginning with the first day of each 
calendar year (except for the year 1947, because the data 
available for this year is from Feb 14 onwards). In this way, 
each time series of one year length can encompass several 
cycles of coronal rotation. The data set is almost continuous 
through out the period of study. We use linear interpolation 
to fill the few data gaps (less than 1%) found in some data 
sets. A few typical examples of such time series are shown 
in the left panel of the Figure [2] The rotational modulation 
due to the solar rotation can be noticed in each of the time 



series. The time series (shown in the left panel of Figure 
[2j are particularly chosen from the years of solar maxima 
(1947, 1969 and 1989) and solar minima (1986, 1996 and 
2009) . The significant rotational modulation of the autocor- 
relation function in most of the years is clear display of few 
nice bright emissive tracers phased well in the longitudes. In 
some cases more pronounced autocorrelation function can 
appear even during minima of the solar cycle, if only one 
isolated (or only very few) active region is on the solar disc 
for several months. 

The autocorrelation function for such time series of the 
solar radio flux at 2.8 GHz of each year is obtained. The right 
panels of Figuref2]show six such plot of autocorrelation func- 
tion out of the 63 obtained in the same manner as discussed 
above. The autocorrelation function shown here corresponds 
to the years of solar maxima (1947, 1969 and 1989) and so- 
lar minima (1986, 1996 and 2009). It can be observed from 
Figure [2] that most of the autocorrelation function curves 
show a very smooth and cyclic rotational modulation with 
a fair amount of correlation. 

The synodic coronal rotation period can be estimated 
from the position of the first secondary maximum. To de- 
termine the synodic rotation period with maximum possible 
accuracy, we fitted a sinusoidal wave whose amplitude is 
the same as first secondary maximum in the autocorrela- 
tion function curve and period accurate to a day. Then we 
considered 7 points (at interval of ±1 day) around the first 
secondary maximum of both the sine curve and autocorre- 
lation curve and calculate the standard deviation. Next, we 
changed the period of sine wave in small steps (in the steps 
of 0.1 day) and tried to minimize the standard deviation. 
The value of wave period for which the standard deviation 
is found to be least would be the synodic rotation period. 
This fitting allowed us to determine the lag with the sub-day 
precision and to estimate the uncertainty of the lag deter- 
mination, which is 0.1 day in present case. 

The synodic rotation periods are then converted into 
sidereal rotation periods using the following expression. 

365.26 x T sv 

i sidereal = \£) 
oDO.^D -t synodic 

This equation does an approximate conversion of syn- 
odic to sidereal rotation period. For more a ccurate con- 
versio n one may use the approach outlined by iRosa et al. I 
(1995). Autocorrelation analysis can thus be used to ad- 
vantage for investigating the periodical variations of time 
series, particularly the long ones. As we know that the in- 
dices of the Sun, derived from the observations or analy- 
sis of various solar activity phenomena, also form time se- 
ries of such type. It is often used with the lag to deter- 
mine the stationarity of a time series. This analysis is useful 
in predicting trends or cyclical functions pre sent in a time 
series. This technique was first employed by lHansen et al.l 
( 1969) to determine the synodic rotation period of the solar 
electron corona after analyzing the K-coronameter obser- 
vations during 1964-1967. Later |E1-Raev fc Scherrer] (|l972T ) 
analyzed sunspot number, solar radio flux and interplane- 
tary field for the period 1967 to 1970 t o det ermine the so- 
lar atmospheric rotation. iParker et all (|l982T ) showed that, 
within the limit of uncertainty, maximum entropy spectral 
analysis and autocor r elatio n analysis provide a lmost similar 
results. iFisher et all l| 19841 ); ISime et al.1 (|l989h applied the 
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Figure 2. A few typical time series of solar radio flux at 2.8 GHz (left panel) for the years of solar maxima (1947, 1969 and 1989) and 
minima (1986, 1996 and 2009) and their respective plots of autocorrelation function used in the present work for the estimation of the 
synodic rotation period (right panel). 
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same analysis technique in their study of the rotational char- 
acteristics of white-light and Fe XIV corona, respectively. 
The differential rotation rates of s oft X-ray features in the 
solar corona were determined by Weber et al. I l| 19991 ) us- 
ing the Lomb-Scargle periodogram. These results were then 
compared with the results obtained through the autocorre- 
lation m e thod and were found to agree within lcr. lVats et al.l 
(199830,1 

200 ll ) used the autocorrelation technique to deter- 
mine the rotation period and the differential rotation as a 
function of height in the solar corona, respectively. 



3.2 Crosscorrelation Analysis 

Crosscorrelation can be used to determine the degree of fit 
between two data sets x and y. The day-to-day measurement 
of annual sunspot number and coronal rotational period, ob- 
tained from the analysis of the 2.8 GHz radio flux of the same 
time interval, can be represented in the form of two time se- 
ries, (350 , 351, 352 , 353, , Xn-l) and (2/0,2/1,2/2,2/3, , 2/n-l)- 

Again, we can form (n — I) pairs of observations from 
these two series. They are: (350,35;), (351,351+;), (x2,X2+i), ■ ■ 
. . ., (x n -i-i,x n -i) and (2/0, 2/i), (2/1,2/1+!)) (2/2,2/2+;): • ■ ■ • 

•, (s/n-/-l,J/n-l). 

The crosscorrelation coefficient P xy {l) of the two time 
series as a function of the lag I is given by, 



Pxy(l) 



Pxy (I) 



ELo 1 ( x k+w -z)(yk -y) 



En — I / 



En — I f 
fc=0 (2/fc 



yf 



ELo 1 (xk-x){y, 



k+i 



V) 



(3) 



(4) 



En — I / 



YJlJoivk-yf 



for I < and I > 0, respectively. Here, x and y is the mean 
of the daily observations of two different time series. 

The cross-correlation analysis helps in determining the 
mutual correlation between solar rotation and other ob- 
served indices of the Sun as well as other stars. Spectro- 
heliograms obtained in extreme ultraviolet (EUV) lines and 
the Lyman continuum are used to determine the rotation 
rate of the solar chromosphere, transition region, and corona 
ijDurjree fe Henzd Il972l ; iHenze fe Dupred Il973n . A cross- 
correlation analysis of the observations indicates the pres- 
ence of differential r otation through the chromos phere and 
the transition region. iBocchialini fe Baudinl (|l995h used this 
method for the diagnostics of the chromospheric dynamics. 
The correct algorithm to calculate the cross-correlation func- 
tio ns of the He I and th e O V line intensities was proposed 
bv lGomorv et all l|2004 ). 



4 RESULTS 

The sidereal rotation period obtained from the time series 
of the radio flux of the years 1947 to 2009 are tabulated in 
Table 13.21 A plot is also given in the Figure (by scattered 
dots), which shows the sidereal rotation period as a func- 
tion of the year. It is found that the mean sidereal rotation 
period during 1947-2009 is equal to 24.3 days. The maxi- 
mum rotation period is found to be 29.5 days in the year 
2007 and the minimum rotation period is 19.0 days in 1970. 



The precision is same in the determination of the sidereal 
rotation period throughout the period of study and is 0.1 
day. The standard deviation of the sidereal rotation period 
during 1947 to 2009 is calculated and is equal to 2.0 days. 

These values of the coronal rotation period show quite 
a large variability (from 19.0 to 29.5 days). Such large 
va riability in the rotation per i od ha s also been reported 
bv lHoward. Gilman fe Gilmanl (|l984h in their study of the 
collection of Mount Wilson white-light plates from 1921 
through 1982. From their Table 2, it can be seen that the 
sidereal rotation rate varies between 11.89 degree/day (in 
1965) to 16.47 degree/day (in 1963) in the equatorial region. 
It means that the sidereal rotation period obtained with 
white-light plates also shows a large variability from 30.3 to 
21.9 days, which is in good agreement with the coronal ro- 
tation period obtained in present work. Ironically, the year 
of maximum and minimum rotation rate (1963 and 1965, 
respectively) are quite nearby, which again shows analogy 
with the present work. In some cases, there is a large vari- 
ation in rotation period in adjacent years (shown in Figure 

GO- 



IChandral l|2010h obtained the coronal rotation period, 
as a function of latitude, using Nobeyama Radioheliogram 
(NoRH) images at 17 GHz for the year 1999 through 2005. 
The sidereal period obtained for the year 2003 shows quite 
a large variability from 29.2 to 22.1 days. This is again 
in agreement with our results obtained analizing radio flux 
data at 2.8 GHz. 

To see the presence of any long term variability in it, 
the coronal rotation period data is smoothen by taking the 
3, 5, 7, 9, . . , 17 years running mean of the sidereal rota- 
tion period. The smoothed rotation periods so obtained are 
then plotted against the year. It is found, after comparing all 
such plots, that the running mean of eleven years smoothed 
the raw data most. By taking the running mean of eleven 
years, the short term variation present in the rotation pe- 
riod data sets are smoothed and thus removed. The data 
obtained after taking the running mean will now show only 
the long term variation, i.e. with the period of more than 
11 years. The plot is shown in Figure [3] (by solid line). An 
evidence of long term periodic variation can be seen in this 
plot. The overall variation has reduced from ~ 11 days to 
only ~ 3 days, as shown in Figure [3] This plot has several 
peaks and dips. However, there are two major maxima and 
two major minima which can be easily identified in Figure 
[3] The separation of consecutive maxima or minima is ~ 22 
years. This is an evidence of ~ 22 years periodicity in the 
temporal variation of sidereal coronal rotation. 

This periodicity may be better seen in the plot of auto- 
correlation function. Hence, we derived the autocorrelation 
coefficient of raw and smoothed (an eleven years running 
mean) of the sidereal rotation period. These are plotted 
against the lag in days in Figure [4] The dashed and solid 
curves are for raw and smoothed sidereal rotation period, 
respectively. The autocorrelation function exhibits the exis- 
tence of a cycle close to a period of 22-years. As expected this 
periodicity is better seen in the autocorrelation of smoothed 
data set. 

The crosscorrelation function between the smoothed ro- 
tation period and sunspots number is calculated and shown 
in Figure [5] by solid line. It is natural that from a crosscorre- 
lation of these time series, we expect to see evidence of any 
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Figure 3. Temporal variation of the sidereal rotation period obtained using autocorrelation analysis of the solar radio emissions at 2.8 
GHz during 1947-2009. The smoothed curve (eleven years running mean) of the sidereal rotation period is also displayed to compare it 
with its raw results, and also to see, if any long term component is present in it. 



Table 1. Sidereal rotation period (in days) obtained through autocorrelation analysis of daily radio 
flux at 2.8 GHz for the years 1947-2009. 



Year 


Sidereal 


Year 


Sidereal 


Year 


Sidereal 


Year 


Sidereal 


Year 


Sidereal 




Rotation 




Rotation 




Rotation 




Rotation 




Rotation 




Period 




Period 




Period 




Period 




Period 


1947 


24.5 
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variation having a period of more than 11 years. The plot of 
crosscorrelation function in Figure [5] (solid line) varies from 
~ —0.30 to 0.53. One can easily identify two maxima and 
two minima in this curve and estimate the period. From the 
positions of these we get a periodicity of ~ 22 years. This 
periodicity in sunspot numbers exists due to the 22-years 
magnetic cycle. The temporal variation of the coronal ro- 
tation period also seems to have a component of the same. 
However, as the central maximum of the crosscorrelation 
function fall very close to lag (at -1 year), the smoothed 



coronal rotation period and sunspot numbers can considered 
to be almost in phase. 



5 DISCUSSION 

Based on the model calculations, we know that this radio 
emission at 2.8 GHz originate from the l ower solar coron a 
near ~ 60, 000 km above the solar surface (|Vats et al.ll200ll ). 
Thus the present statistical analysis of solar flux at 2.8 GHz 
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Lag, in year 

Figure 4. Autocorrelation function curves of raw and smoothed 
annual sidereal rotation period for the period 1947-2008. 




Lag, in year 

Figure 5. Crosscorrelation function of smoothed data of annual 
sidereal rotation period and sunspots numbers. 

reveals some new long term rotational features of the so- 
lar corona at a height ~ 60, 000 km above the photosphere. 
The results obtained through the autocorrelation and cross- 
correlation analysis can be discussed in the following three 
parts. 

(i) The temporal variation of the rotation period (Figure 
O does not show any systematic periodicity. On the other 
hand sunspots have a very clear 11 years periodicity. So 
from this, it seems that probably there is either no or very 



weak correlation between the coronal rotation period and 
the sunspot cycle. 

(ii) A reasonable presence of a long term variation in the 
coronal rotation period ( of ~ 22-years) become visible sim- 
ply by smoothing (11 point running mean) the values (Fig- 
ure[3]). This may correspond to the 22-years magnetic rever- 
sal cycle or Hale's cycle. 

(iii) The crosscorrelation analysis between the smoothed 
data of coronal rotation (Figure [3} and its correlation with 
the sunspots number also shows the clear presence of a 22- 
years Sun's magnetic field reversal cycle (Figures [4] and [5} . 

The position of minima in the Figure reveals that the 
years of minima of the 22-years cycle, obtained after smooth- 
ing the rotation period data, coincide well with the years of 
minimum activity of even numbered solar cycles (ENSCs), 
i.e. at the years 1954, 1976 and 1996 (shown in Figure [TJ. 
Whereas, the years of maxima of the 22-years cycle matches 
well with the years of minimum activity (at the years 1964 
and 1986) of odd numbered solar cycles (ONSCs ). This re- 
sult is in agreement with the results obtained by I Javaraiahl 
(2000). The variation in equatorial rotation rate A is found 
to be significant only in the years of ONSCs. Where A is re- 
ciprocally related to the equatorial sidereal rotation period 
(T) by an expression (T = 360 /A). It can be seen again in 
Figure [3] (solid line) that the smoothed rotation period rises 
steeply in the ONSCs and falls back gradually in the ENSCs. 
A similar trend continues even in solar cycle 23. This shows 
a good agreement in both the results obtained through two 

entirely different methods and tracers. 

The Magnetic butterfly diagram due to lHathawavl 

(1998) shows clearly that the reversal of the Sun's mag- 
netic polarity is triggered in the years close to the years of 
solar activity maxima, i.e. years 1979, 1989 and 2000. But 
Figure [3] indicates that the ascent or descent in the value 
of the smoothed sidereal rotation period is initiated in the 
years of solar minima. Therefore, it shows that there may be 
anticorrelation between the magnetic pole reversal and the 
rotation period. 
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